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Abstract 

The role of collective modes in various experiments on the cuprates is in- 
vestigated. We calculate the neutron scattering, photoemission (ARPES), 
and Raman scattering intensities below Tc within the fluctuation-exchange 
(FLEX) approximation for the two-dimensional Hubbard model. It is shown 
that the large peak in the dynamical spin susceptibility arises from a weakly 
damped spin-density-wave collective mode. This gives rise to a dip between 
the sharp low energy peak and the higher binding energy hump in the ARPES 
spectrum. Furthermore, we show that the collective mode of the amplitude 
fluctuation of the d-wave gap yields a broad peak above the pair-breaking 

threshold in the Big Raman spectrum. 
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A wide variety of collective modes has been observed in the three phases of superfluid 
These fluctuations of the spin-triplet p-wave gap have been calculated from cou- 
pled Bethe-Salpeter equations for the T-matrices in the particle-particle and particle-hole 
channels0i. This method has also been used to investigate the collective modes in hypothet- 
ical p-wave pairing superconductors like Heavy Fermion sup er conduct or s@. A detailed study 
of the collective modes for 3D d-wave superconductors, including different order parameter 
symmetries, has been made in Ref. |. Recently, pair fluctuations and the associated Raman 
scattering intensity have been calculated for a two-dimensional (2D) li-wave weak-coupling 
superconductoiil. 

In this note we investigate the collective modes within the fluctuation-exchange (FLEX) 
approximation for the 2D one-band Hubbard model and their relevance to neutron scatter- 
ing, photoemission, and Raman scattering experiments in high-T^ cuprates. The FLEX- 
approximation for the particle-hole channefl^ yields the dynamical spin susceptibility, 
Xs = Xso(l - Uxsoy^, and the charge susceptibility, Xc = Xco(l + Uxco)~^- Here, U is 
the on-site Coulomb repulsion, and Xso{fi,^^) and Xco{fi,^) are the irreducible susceptibili- 
ties. The latter are calculated from the dressed normal and anomalous Green's functions G 
and F. The corresponding normal self-energies, uj[l — Z(]s.,uj)] and ^(k, cj), and the d-wave 
gap function, 0(k, cj), are determined self-consistently by the Eliashberg equations with in- 
teractions given by (3/2)f/^Im Xs{^.,^) =t (l/2)t/^Im Xc{^,^) (plus (minus) sign for normal 
(anomalous) self-energy). 

Below Tc large peaks evolve in the spectral density Im Xs{<i,^), i-e., four distinct peaks 
at wave vectors q near Q = (tt, vr) for next-nearest neighbor hopping t' = oi, and a broad 
peak centered at Q for t' = — 0.45t (t is the near neighbor hopping energy )§. These results 
are in qualitative agreement with neutron scattering experiments on La2_xSr^Cu04 and 
YBa2Cu307_5, respectively^. Similar results have been obtained within the t-J mo del0. 
In Fig. |I] we show Im Xs(Q,^^) for U = 3.6t, t' = — 0.45t, and band flUing n = 0.90. One 
sees that a large peak evolves at about ujq = 0.08t as T decreases below = 0.022i(:. The 
amplitude Aq of the (i^2_y2-wave gap rises much more rapidly below Tc than the BCS d- 



wave gap and reaches at our lowest temperature T = 0.017t {T/Tc = 0.77) a value of about 
Aq = O.lt (see Fig. |^). We find that the peak in Fig. |l]is due to a slightly damped collective 
mode because the susceptibilty has a pole at c^o, more exactly, Re Xso{Q,^q) — U^^ = 0, 
and the height of the peak is large of the order of the quasiparticle lifetime l/r(ci;o)- Here, 
r(k, cu) = culm Z(k, C(j)/Re Z(k, cj) is the quasiparticle scattering rate. Since this is decisive 
for the observability of the collective modes in the cuprates we show in Fig. ^ the functions 
ulm Z(k, cij) and Re Z(k, cj) at the anti-node k^ and the node k^ of the gap on the Fermi 
line. One sees from Fig. |^ that for T below Tc the scattering rate decreases dramatically for 
frequencies u below the pair-breaking threshold 2Aq ~ 0.2t. 

In order to understand somewhat better the origin of this spin-density-wave collective 
mode we have calculated Xso(Q)^) in the weak-coupling limit. The sums over Matsubara 
frequencies have been carried out with the help of the methods developed for superfiuid 
The effect of quasiparticle damping is taken into account by carrying out the analytical 
continuation of this result from to u + iT. For a gap A(k) = (Ao/2)(cos A;^ — cosA;^^) 
and a band e(k) with t' = and chemical potential /i the summation over k in the square 
Brillouin zone has been carried out numerically in the following expression for T = 0: 

■Sp EkEk+Q — (-k^k+Q — ^k^k+Q Ek + Ek+Q 

V {Ek + Ek^Qf - + iTf 2EkEk+Q ■ ^' 
Here, = e^(k) + A^(k). Then we obtain a peak in the t^-function Re XsoiQ,^) at the 
kinematical gap u = 2|/i|Ei whose height decreases with increasing F. The approximate 
analytic result for T = is given by 

X.o(Q,^) = V,-' - Np{z/1 + z)'/Hog [4 (1 + zf'] ; 

z = [V - + ^r)'] / (2Ao)' ; 
yo-' = iV^log(2iy/Ao) . (2) 

Here, W = At is the half bandwidth. The function Re Xso{Q,^^) in Eq. (|^) rises first 
with u;^ and then exhibits a peak at the kinematical gap 2\fi\ whose height is about V"o~^ — 



(7r/2)A'"i?(r/2Ao). A low-frequency mode, i.e., a zero of equation Re XsoiQ,^) = '^/U, is 
obtained only for a finite range of U values which decreases with increasing F. For t' = — 0.45t 
a kinematical gap no longer exists and the effective |/i| is nearly zero. Then the approximate 
analytical result for the expression in Eq. (0) becomes equal to: 

1 

XsoiQ, ^) = V^o"' + -Npi———K{u + ij) ; 
2 [uj + 27) 

cD = cu/2Ao ; 7 = r/2Ao . (3) 

Here, K is the first elliptic integral. Now the peak of Re Xso as a function of uj occurs 
approximately at ~ 2Ao. This result has been checked by carrying out numerically the 
sum over k in Eq. ([l|) for t' = — 0.45t and for different amplitudes Aq and chemical potentials 
/X. In fact, we find that the function Re XsoiQ,^) exhibits a peak at about u ~ 2Ao 
whose height decreases as fi increases, for example, from —1.3 (unrenomalized band filling 
uq = 0.84) to —0.8 {tiq = 1.03). A solution uq < 2Aq of the equation Re Xso(Q)'^o) = 
exists only for a small range of U values near U ^ 3t. The strong-coupling calculation 
yields a smaller resonance energy. Another difference in comparison to the weak-coupling 
result is the fact that the self-consistent strong-coupling calculations yield a collective mode 
for much higher values of U, for example, U = 6.8t in Ref. This shows how important 
it is to take into account the feed-back effect of the self-energy on the dynamical spin 
susceptibility Xs- These strong renormalization effects might also be responsible for the 
observed broadening and decrease of the resonance energy of the neutron scattering peak in 
underdoped YBa2Cu306+x' which goes in proportion to the decrease of Tc or doping level0. 
In fact, for decreasing doping x = n — 1, or increassing chemical potential /i, the position 
and height of the function Re Xso(Q) ^) decrease which means that the position of the peak 
of Im Xs is decreased and its width is increased. 

We show now that the spin- density- wave collective mode has a large effect on the angle- 
resolved photoemission intensity (ARPES) below Tc. In Fig. ^(a) we have plotted our results 
for N{'k,uj)f{uj) (where A^(k, cu) is the quasiparticle spectral function and f{uj) is the Fermi 
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function) for several k-vectors ranging from k = (vr, 0), (77r/8,0), (137r/16,0), . . ., down to 
(0, 0). The parameter values are the same as in Figs. |I] and |^. For k-vectors near (tt, 0) we 
have a sharp low energy peak followed by a dip and then a hump at higher energy. For k 
moving from (tt, 0) to (0, 0), the sharp peak remains first at about the same position while the 
broad hump moves to higher binding energy. In Fig. ^(b) we show the corresponding normal 
state spectra at T = 0.023t. One notices that the broad hump at higher binding energy 
remains at the same position upon entering the normal state, while the sharp peak and the 
dip feature disappear. In the superconducting state along the nodal direction of the d-wave 
order parameter we do not observe the dip feature as we have shown previously^. These 
results are in qualitative agreement with the photoemission spectra of Bi2Sr2CaCu308+5 (Bi 
2212)0. In this paper it was argued that the dip in the spectrum stems from a step-like 
edge in the quasiparticle scattering rate which arises from the interaction with a collective 
mode. This scenario is confirmed by our results for the collective mode shown in Fig. |I| and 
by the scattering rate shown in Fig. ^ We estimate from the edge of the peak in Fig. |^ 
a gap amplitude Aq — O.lt at T/Tf. = 0.77 and a spectral dip at binding energy of about 
2.3Ao corresponding to a mode frequency Uq ~ I.SAq according to the estimates of Ref. p . 
However, here we have a discrepancy with regard to the latter estimates because our mode 
frequency shown in Fig. [l| is much lower, i.e., uq ~ O.SAq. We note that we obtain also 
a dip in the density of states N{uj) below the gap peak at negative uj values which agrees 
qualitatively with STM measurements on Bi 22120. 

We want to mention that higher order peaks in the photoemission spectra due to the 
collective mode as have been observed for example in solid hydrogen0 are not visible here. 
This is due to the fact that in our case the spin-density-wave collective mode is a damped 
mode and the strong quasi-particle damping rate washes out higher order peaks. In addition, 
the self-energy contains an average over momentum, further reducing this effect. This is 
documented by the fact that the normal state spectrum in Fig. ^b), where there is no 
collective mode present, is not much different from the spectrum in the superconducting 
state in Fig. §(a) at higher binding energy. 



We come now to the discussion of order parameter collective modes in (i-wave super- 
conductors which can be calculated in analogy to those in p-wave pairing superconductors!. 
In general it can be said that the (i2.2_j^2-wave pairing component in weak-coupling theory 
gives rise to the phase fluctuation mode which is renormalized into a 2D plasmoni, and 
to the amplitude fluctuation mode of the (i-wave gap. For each additional (weaker) pair- 
ing component, like an extended s-wave component, one obtains an amplitude (real) and a 
phase (imaginary) fluctuation mode. Let us flrst consider the amplitude fluctuation mode 
of the dx2-y2-wave gap. We have calculated the mode frequency uq from the weak-coupling 
expression in Ref. for q = 0: 



Re 



X 



^(cj2 - 4A^) [cos(A;^) - cos{ky)] 

L k 

1 tanh(Efc/2T)" 



AEt -iu; + try 



(4) 



By summing numerically over k in the square Brillouin zone we obtain for t' = in e(k) 
and T = two solutions with frequencies ujq ~ a/3Ao provided that the damping V is 
sufficiently large, namely, uo^ < 3.5r. For t' = — 0.45t and T = we obtain two solutions 
whose frequencies are somewhat larger, ujq ~ 2Ao, where again the condition ujq < S.SF has 
to be satisfied. For a mode frequency = 2Ao — 0.2t at T/Tc = 0.77 (see Fig. H) one finds 
from Fig. ^ a damping r(ti;o) — O.lt at the anti-node k^ which means that the condition 
ujq < 3.5r is satisfied. In Ref. ^ a frequency ujq = VSAq was obtained for the amplitude 
collective mode, however, the coupling of this mode to the charge fiuctuations was neglected. 
We find that the coupling of this fiuctuation in the particle-particle channel to the charge 
fiuctuation in the particle- hole channel yields approximately the following contribution Xfi 
to the charge susceptibility Xco at T = (see Refs. ^ and ^): 



(5) 



2-2 4 2 , 8._ 

-uj H — 7 — 1 tuj'y 

3 3 ' 3 ' 



+ 7 (Aco'^ - 27^ + Gito'j^ log (4 1 - (tu + 27) 



-1/2 



cD = cu/2Ao ; 7 = r/2Ao . (6) 

Here, Np and Np = dNp / du are the density of states and its derivative at the Fermi energy 
C(j = 0. One notices from Eq. (BT) that in the hmit 7^0 one obtains no vahd solution 



of the equation Re g{ujQ) = because the solution ujq = y3/2 violates the condition that 
cJo < 1- However, for suffiently large values of F (7 > a;/2) one obtains a solution of the 
equation Re g{uJo) = which satisfies the condition uJq < 1. One sees from Fig. || that 
this condition is approximately satisfied for ujq = 2Ao — 0.2t because then one enters the 
pair-breaking continuum where F ~ co'/2 near the anti-node k^. Thus the mode frequency 
is about UJq = 2Aq for damping F = Aq in agreement with the numerical results. 

We have calculated the resonance frequency of the exciton-like s-wave mode of the order 
parameter which is caused by an additional s-wave pairing component l^fol which is smaller 
than the main d-wave pairing component \g2\ (see Ref. The method of Refs. ^ and ^ 
yields the following contribution Xexc of this order parameter fluctuation mode to the charge 
susceptibility Xco at T = 0: 



Xe 



0,u;) = -{Npujf[g,,,{uj)r' , (7) 



90/ k '^^k 
1 tanh(igfc/2r) 



2 V Ek [4El -{u + tT 
From Eq. (H) we obtain the following approximate result: 

9eU^) = f 1 - -1 ^ + Npt—^K{LO + 27) (9) 

where V^~^ is given by Eq. (Q) and ij and 7 by Eq. (^). We have carried out the summation 
over k in Eq. (^) numerically and find in agreement with Ref. ^ that a solution of the 
equation in Re gexci^o) = at given Aq and F exists only for very small values of the 
parameter ((72/fi'o) ^ 1 0.1). This means that the s-wave pairing coupling has to be 
almost as strong as the rf-wave pairing component which is quite unrealistic. However, for 



increasing T the resonance frequency ujq decreases and becomes much smaller than the pair- 
breaking threshold 2Ao for reasonably large scattering rates (r/2Ao ~ 1/2)- This means 
that the contribution Im Xexc{^) of the exciton-like mode to Raman scattering intensity 
with Big polarization shows up as a small peak below the pair-breaking threshold. Since 
the damping F in the direction of the momentum of the antinode of the order parameter 
rises rapidly with u (see Fig. §(a)) it may be that this peak becomes observable for smaller 
values of the ratio qq / §2 of s-wave and (i-wave pairing couplings than those which have been 
obtained from weak-coupling theory!. 

In the weak-coupling limit it has been shown that vertex corrections due to the d-wave 
pairing interaction together with electron-electron scattering lead to good agreement with 
the Big Raman data on YBCO0. In Fig. HI we show our strong-coupling results for the 
Raman response functions Im x-yisi = 0) ^) where 7 are the vertices 7 = t[cos{kx) — cos{ky)] 
and 7 = —At' sm{krc) sin{ky) for Big and i?2g symmetry. One sees that for Big symmetry a gap 
and a pair-breaking threshold develop below with a threshold at about 0.15t (3/2)Ao 
at T/Tf. = 0.77 (see Fig. |^). This means that the peak of the order parameter collective 
mode at ujq ^ 2Ao and width Aq lies in the pair-breaking continuum. The question arises 
whether or not the contribution of Im Xfi to the Big Raman spectrum is sizeable because the 
coupling strength proportional to N'p/Np in Eq. (^) arising from particle- hole asymmetry 
is rather small. However, in the strong-coupling calculation the coupling strength of this 
mode to charge density given by TY,k^nG(}^i'>'^n+m)F{\<i,iuJn) is much larger. The reason 
is that beside the term proportional to e(k) yielding N'p/Np, one obtains additional terms 
proportional to the self-energy components Re ^(k, cu) and Im ^(k, cu) which give relatively 
large contributions. In addition, one obtains a contribution from the imaginary part of the 
gap function, i.e. Im 0(k, cj). 

In conclusion we can say that the spin-density-wave collective mode below gives rise 
to large effects in the magnetic neutron scattering and photoemission intensities and the 
tunneling density of states. In order to explain the physical basis of our strong-coupling 
results we have compared them with analytical expressions derived from weak-coupling 
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theory. This shows that the gap in the scattering rate and the strong mass enhancement 
of the quasiparticles below are decisive for the observabihty of this mode. On the other 
hand, the amphtude fluctuation mode of the d-wave gap couples only weakly to the charge 
fluctuations and yields a broad peak above the pair-breaking threshold in the Big Raman 
spectrum. This peak may be, at least partially, responsible for the observed broadening 
above the pair-breaking peak because the coupling strength due to particle-hole asymmetry 
is enhanced by strong-coupling self-energy effects. 

Previous work on collective modes in high-Tc superconductor^l&ill has been restricted 
to weak-coupling and mean-field calculations. The FLEX approach we use here, is a self- 
consistent and conserving approximation scheme, which goes well beyond mean-field. Es- 
pecially the feed-back effect of the one-particle properties on the collective modes in the 
superconducting state is included self-consistently and the importance of the quasiparticle 
damping becomes clear. It is therefore a highly non-trivial and satisfactory result, that the 
resonance in the spin-susceptibility, the step-like edge in the quasiparticle scattering rate, 
and the dip-features in the ARPES and tunneling spectra can all be understood within one 
theory in a self-consistent fashion. The self-consistent calculation also yields a larger cou- 
pling strength of the d-wave amplitude mode to the charge density and a lower resonance 
frequency of the s-wave exciton-like mode of the order parameter which makes it more likely 
that these modes show up in the Big Raman scattering channel. 
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FIGURES 

FIG. 1. Spectral density of spin susceptibility at wave-vector Q = (tTjTt), ImXs(Q,a;), for 
temperatures T = 0.023t, 0.020i, and 0.017t (Tc = 0.022t). Here, U = 3.6t is the on-site Coulomb 
repulsion, t the near neighbor hopping energy, t' = —0.45^ the next-nearest neighbor hopping, and 
n = 0.90 the renormalized band filling. 

FIG. 2. Quasiparticle scattering rate, r(k, a;) = ulm Z(k, a;)/Re Z(k, a;), at anti-node and 
node kf, of the d-wave gap, in the normal state at T = 0.023t (solid lines), and in the supercon- 
ducting state at r = 0.017i {T/Tc = 0.77) (dashed lines), (a) ulm Z'(k, w); (b) mass enhancement 
Re Z{k,Lo). 

FIG. 3. Photoemission intensity, N (k, uj) f (u) (here N is the quasiparticle spectral function 
and / the Fermi function) for k = (A;7r,0) where k = 1, 7/8, 13/16, 3/4. 5/8, 1/2, 3/8, 1/4, 1/8, 
and 0. (a) in the superconducting state at T/Tc = 0.77. The narrow peaks at low binding energy 
decrease and vanish, and the binding energies of the broad humps increase in the sequence of k 
values, (b) in the normal state at T = 0.023t. Note, that the broad humps are at the same position 
as in the superconducting state. 

FIG. 4. Raman spectra Im x-yi'i = 0)<^) for Big symmetry at T = 0.023i (solid line) and 
T/Tc = 0-77 (dashed line), and Raman spectrum for B2g symmetry at T/Tc = 0-77 (dotted line). 
Tc = 0.022i. 
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